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Abstract 

We study the quark mass function on hypercubic lattices, in a large 
range of physical volumes and cutoffs. To avoid the very large Wilson 
term artefact, we exploit the relation between the quark mass function 
and the pseudoscalar vertex in the continuum. We extrapolate to the 
chiral limit. 

In function of the physical volume, we observe a striking discontinuity 
in the properties of chiral extrapolation around a physical volume L c ~ 
6 GeV -1 — 1.2 fm. It is present in the quark mass function, which 
collapses to zero, as well as in the pion mass and the quark condensate 
as directly calculated from the pseudoscalar correlator. It is strongly 
reminiscent of the phenomenon of chiral symmetry restoration observed by 
Neuberger and Narayanan at Nc = oo around the same physical length. 

In the case of spontaneous symmetry breaking, we confirm that the 
OPE of the quark mass function, involving the quark condensate, is not 
operative at the available momenta, even taking into account the unusu- 
ally large high order corrections to the Wilson coefficient calculated by 
Chetyrkin and Maier ; the gap remains large, around a factor 2, even at 
the largest momenta available to us (p ~ 6 GeV). 



1 Introduction 

Obviously, when chiral symmetry is spontaneously broken, the quark pseu- 
doscalar vertex presents a Goldstone pole, i.e. a pole at q 2 = m%. In the 
Euclidean region tested by lattice QCD, the large magnitude of its contribution 
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to the pseudoscalar vertex at q = was noted first in numerical calculation by 
ref. pp. However, the denominator is q 2 + rn 2 ; therefore it vanishes at q = 
in the chiral limit and it is proportional to the current quark mass ; we have 
then by exception a divergence of the vertex in the region tested by lattice 
QCD, i.e. it goes to infinity when m q —> 0, for all momenta p of the quark legs 
(q = p' — p = 0), a phenomenon which would be very spectacular if we could 
indeed approach very small quark masses on the lattice. One must note that in 
spite of this divergence, a chiral limit of the vertex was presented in the litera- 
ture for quite a while in the context of MOM non perturbative renormalisation 
at zero mass 2 . Attention to this divergence was drawn by the JLQCD collabo- 
ration 2 , and around the same time in [3l|4] . In the paper [3] , the problem was 
shown to have two distinct aspects : 

1) Usually, the pseudoscalar vertex is used to renormalise the pseudoscalar 
density in the UV asymptotic region ; one uses first a numerical MOM renor- 
malisation constant Zp IOM to be later inserted into certain perturbative calcu- 
lations (running to higher scales, scheme conversions). It is then necessary to 
subtract the Goldstone boson pole, in order to extract the purely perturba- 
tive part of the renormalisation constant. Analogously, one should extract also 
any power corrections originating in condensates, which will be the case for the 
quark condensate contribution, when determining Z m , see below, or for the A 2 
condensate contribution in the gluon and ghost Green functions when extract- 
ing perturbative informations like Aqcd, see for instance |28j and references 
thereins. 

2) However, as a second aspect, one may also consider the Goldstone pole 
contribution not as a parasitic contribution to be eliminated, but in its own 
physical interest and in its relation with other interesting quantities. We recall 
these relations, stressed in the same paper. First, the residue of the pole is 
simply given, up to a finite renormalisation, by the quark mass function, i.e. 
the quotient of the scalar part of the quark propagator and the vector part, in the 
chiral limit. That this mass function is non zero in the chiral limit is therefore 
also a signal of spontaneous symmetry breaking. This can be seen directly from 
the fact that the scalar part of the propagator is zero in a chirally symmetric 
vacuum. Then, this mass function can be said to constitute a continuous set 
of renormalised order parameters, labelled by p. 

Finally, the spontaneous breaking also manifests itself through the OPE of 
this quark mass function ; the main contribution in the chiral limit is the one of 
the quark condensate, with power 1/p 2 ; this gives therefore a relation with 
another, well known, indicator of spontaneous chiral symmetry breaking. 

The object of the present paper is to study in more detail these relations 
on the lattice, at Nf = 0, with the clover action. The central object will be 
the quark mass function, in the chiral limit, calculated, along the line described 
above, through the residue of the Goldstone boson in the pseudoscalar ver- 
tex. Thus the present paper completes the study of the papers OHO] devoted to 
the vector part of the propagator, For an extensive study of the quark prop- 
agator, one can refer also to the papers of the Adelaide group, see for instance 
our reference below [IT] . 

Although it may appear involved at first, this method of calculation is indeed 

2 Let us stress that no such Goldstone contribution is expected in the method of renormal- 
isation of the ALPHA group at zero mass. 
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an advantageous approach in calculations with the clover action, because the 
lattice quark propagator itself would yield an untractable mass function, plagued 
with the huge Wilson term artefact, while the pseudoscalar vertex has much less 
artefacts, as will be seen. 

On the other hand, this calculation is very simple and fast even on large 
volumes in contrast to the calculations with overlap action. Of course, one 
cannot claim to work very close to the chiral limit. It is however instructive to 
see what comes out when using the standard methods of chiral extrapolation, 
successfully used since a long time in phenomenological applications of lattice 
QCD. 

With respect to the exploratory study presented in Ref. ( [3l|4]) with the 
help of data from the QCDSF group, the present study uses the same basic 
ideas, but it is performed with better data, with various /3's and lattice sizes, so 
that for instance artefacts can be identified and OPE can be tested with larger 
momenta. 

1.1 Physics results 

OPE : let us recall, with the lattice data used in [3], there appears a very 
large discrepancy between the prediction of OPE and the lattice data, using the 
known quark condensate value and the Wilson coefficient with the two calcu- 
lated orders ; the prediction of the OPE is found to be much lower than the 
lattice data. This problem with OPE is confirmed (see section [5]) • It can be 
partly elucidated thanks to the new perturbative calculations of Chetyrkin and 
Maier for the Wilson coefficient [5j[6] : the high order corrections to the Wilson 
coefficient are exceedingly large, almost spoiling the hope to extract the con- 
densate from the high momentum lattice data. The discrepancy is admittedly 
reduced, but remains around a factor 2 at the largest available momenta. 

Chiral symmetry "restoration" : one also displays (see section 0]) a 
new striking and unexpected phenomenon at small volumes, that we can term 
for simplicity as an abrupt "chiral symmetry restoration" below some critical 
physical volume ; it affects the chiral extrapolation of several quantities at the 
same time, not only the quark mass function, but also the pion mass and the 
quark condensate ; to evaluate the latter independently of the GMOR (Gell- 
Mann, Oakes, Renner) identity, we recourse to the old method proposed in the 
paper by Bochicchio et al. [7] . This discontinuity is reminiscent of the phase 
transition found by Narayanan and Neuberger [5] at N c — > oo. Let us recall at 
this point that usually it is expected that the chiral symmetry order parameters 
are vanishing at zero quark mass, in finite volumes ; therefore, what one should 
observe on the lattice is depending on the physical volume and on the smallest 
quark masses actually considered. 
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2 Setting the framework in the continuum : the 
pseudoscalar vertex, the W-T identity and the 
quark mass function 

2.1 Definitions and Lorentz invariance 

We work in the Landau gauge. Let us first fix the notations that we will use. 
We will use all along the Euclidean metrics. The continuum quark propagator 
is a 12 x 12 matrix S(p) for 3-color and 4-spinor indices. One can take into 
account Lorentz (in fact 0(4)) invariance and discrete symmetries, as well as 
color neutrality of the vacuum by expanding the inverse propagator according 
to : 

S- 1 { P ) = 5 a , b Z^{p 2 )(ij> + m{p 2 )) (1) 

where a, b are the color indices. Z^{p 2 ) is a standard lattice notation , referring 
to the role it plays as a renormalisation constant for the quark field (in the stan- 
dard Georgi-Politzer MOM renormalisation, see subsection 12.31 Z^ 1 OM (/z 2 ) = 
Z-ipdi 2 ) ; for the precise lattice definition, sec below, section [3]). Obviously, one 
has in the continuum, with trace on spin and color : 

Z^(p a ) = l/il/12Tr(S- 1 (p) J J)/p 2 (2) 

Sometimes, one uses to describe the scalar part of the propagator, instead of the 
quark mass function m(p 2 ), the alternative quantity, which is the proper scalar 
part of the bare propagator : 

b{p 2 )^Z^{p 2 )m{p 2 ). (3) 

One must not forget however the advantage of m(p 2 ) being UV finite, since 
it is the ratio of two quantities renormalised by the same factor Z^. In fact it 
is the MOM renormalised quark mass, with renormalisation point /j? = p 2 (see 
below) . 

This is the central physical quantity we will consider in the present paper. 
On the other hand,we will consider as auxiliary the pseudoscalar vertex, and 
then use the axial Ward identity to relate this vertex to m(p 2 ). Indeed, it 
happens that m(p 2 ) is not directly calculable with the clover action. We then 
need some more definitions. 

Let us consider a colorless local two quark operator qOq. The corresponding 
three point Green function G is defined by 

G(p, q)= [ d*xd*y e^y+^ x < q(y)q(x)Oq(x)q(0) > (4) 

It is a 4 x 4 matrix in Dirac space. The corresponding vertex function is then 
defined by amputation of quark propagators on both sides : 

r(p,q) = S- 1 (p)G(p,q)S- 1 (p + q) (5) 

In the whole paper, we will restrict ourselves to the case where the operator 
carries a vanishing momentum transfer q^ , = 0. In the following we will omit to 
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write = and we will moreover understand as the bare vertex function 
computed on the lattice. 

Now, Lorentz covariance and discrete symmetries allow to write for the axial 
vertex : 



9a\p 2 )p^15 + *.9a 4) (p 2 )[7ao$75] 



(6) 



which should be obeyed approximately on the lattice, as we checked, and simi- 
larly for the pseudoscalar vertex :. 



r { P z ) = S a , b g < i\p 2 )l 5 + ^g^(p ^ )l^ 



(7) 



2.2 Re normalisation and Ward-Takahashi Identities 

Although we do not require a specific renormalisation scheme, we have to dis- 
cuss the renormalisation, because the Ward-Takahashi(W-T) identities should 
be imposed on the renormalised theory, and not on the bare quantities (wc 
do not consider anomalies). The following considerations hold in an arbitrary 
renormalisation scheme, and the corresponding renormalised quantities are de- 
noted by a sub- or superindex R. The renormalised equations should hold up 
to 0{a) artefacts. 

Z 2 denotes as usual the fermion field or propagator renormalisation according 

to : 

q = \/Z^qn 

S( P ) = Z 2 S R ( P ) (8) 
Let us recall that the corresponding renormalised vertex functions arc : 

r(p) = z 2 - 1 z 1 r R ( P ), (9) 

where the necessary subindices are implicit for each type of vertex ; Zq is 
the renormalisation of the composite operator, namely a current or density 
operator : O — jv,jA, P5 ; the Z 2 factor takes into account the amputation . 

Note that the standard definition of renormalisation constants is to divide 
the bare quantity by the renormalisation constant to obtain the renormalised 
quantity (except for photon or gluon vertex renormalisation factors Zi which 
we do not use). In principle, renormalisation of composite operators, for in- 
stance Zy, should be defined similarly. We have followed this convention in our 
works on gluon fields, for the renormalisation of the gauge dependent gluon field 
condensate A 2 . But, in the case of quark composite operators, an opposite con- 
vention has become standard in lattice calculations : (qOq)bare = Z^{qOq)Yi \ 
we feel compelled to maintain this convention for the sake of comparison with 
parallel works on the lattice. This explains our writing of the renormalised 
vertex function. 

In the continuum limit Zy = 1, and in the chiral limit, Za = 1 (conserved 
currents). We keep Za since the axial current is not conserved away from the 
chiral limit, which we take only in the end. The lattice artefacts, as any other 
rcgularisation scheme, generate finite 0(g 2 ) effects, vanishing slowly with a, 
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included in factors Zy, Za, due to additional divergencies multiplying the a 
terms (which have higher dimension). There are also terms with powers of a 
which we do not write. Zy, Za are independent of the renormalisation scheme 
up to such terms. The fact that we do not include such terms means that our 
equations should hold only sufficient close to the continuum. 

Consequences of the axial Ward identity Let us develop the conse- 
quences of the axial W-T identity, which derives from the equation <9 m (ja) p = 
2mP§. For this purpose, one has to return momently to the general case 
q = p' — p 7^ 0. Since they reflect the symmetries of the physical theory, the 
naive Ward identities should a priori hold for the renormalised Green functions 
(except for anomalies) and at infinite cutoff, which means : 

<? M I^rCp, g) = -i(5 R ; 1 (p + (7)75 +75<S' R ; 1 (p)) +i2mnT 5 ,n(p,q) (10) 

tor is the renormalised mass in the scheme, to-r = Z~^vn q . It is then possible 
to return to bare quantities which are the ones actually measured on the lattice. 
Multiplying both sides by Z^ 1 one gets : 

Z A «T^(p,flf) = -i(S- 1 (p + q) l5 + l5 S- 1 (p))+i2m R Z P T 5 (p,q), . (11) 

i.e. there appear the renormalisation constants Za,Zp. Now, we exploit this 
equation through an expansion in powers of q^. 

Setting q = 0. The pseudoscalar vertex. First, from (HJ)-©, one finds 
constraints for the pseudoscalar vertex by making q^ = 0, which eliminates the 
axial vertex : 

m(p 2 )Z^,(p 2 )-f 5 = m R Z P T 5 (p) (12) 

This means that the pseudoscalar vertex contains only the 75 component, i.e. 
the second term in ([7]) should vanish : 

<?fV) = o 

r 5 (p) = ffi 1) (p 2 )7 5 . (is) 

This vanishing of gi 2 ^ (p 2 ) is verified to a good accuracy on the lattice. Moreover, 
the pseudoscalar vertex at q — is entirely determined from the mass function 
and Z^,(p 2 ), i. e. from the propagator, through eq. (|T2|) . if we know the propor- 
tionality constant torZp. However, this is not a practical way to determine the 
pseudoscalar vertex on the lattice, because of the very large Wilson artefact in 
m(p 2 ). Rather, as we propose, this relation should be used in the reverse way : 
to determine m{p 2 ) from the pseudoscalar vertex and Z^{p 2 ). 

The proportionality constant torZp is obtained in a familiar way ; we de- 
fine a bare mass p through the equation <9 M j^ = 2pP$, which we can easily 
measure on the lattice through bare matrix elements involving these opera- 
tors (in the technical lattice practice, we rather denote by p the same quan- 
tity in lattice units). Comparing this equation with the renormalised one 
dfj, {ja ) r = 2to r (P 5 ) r , we get at once : 

torZ p = Z A p- (14) 

The good point in this transformation is that the r.h.s is independent of the 
renormalisation scheme, because it is so for Za as we will show now. 
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Indeed, from eqs. (fT2|) and ([LP) , one gets : 

m {p 2 ) = Z A pg ( i\p 2 )/Z^) (15) 

Note that this relation has been demonstrated independently of any specific 
renormalisation scheme. Then Za can be expressed as the same definite com- 
bination of bare quantities Za = b{p 2 ) / p / g\p{p 2 ) in all schemes. However, 
this expression is not a practical way to determine Z A on the lattice, once again 
because of the very large Wilson artefact in the quantity b(p 2 ) defined in eq. ([3]). 

What is now to be noticed is that neither side of eq. f|l 5f) vanishes in 
the chiral limit if there is spontaneous breaking of the symmetry. As 
to the m(j> 2 ) side, this is because the vacuum is then not invariant, therefore 
the expectation value of a scalar need not vanish, while on the other side there 
is a Goldstone pole in g^\p 2 ) compensating for the vanishing of the p factor. 
We then write : 

m chira i(p 2 ) = lim Z A (m q )p g^(p 2 )/Z^(p 2 ) = 
m— s-0 

Z A (m q =0) hm pg^ip^/Z^p 2 ) (16) 

This is the basic equation which we use below to deduce m(p 2 ) in the chiral limit. 
We need both the functions Z^p 2 ) and g^(p 2 ), which we determine from the 
measurement of the propagator and pseudoscalar vertex, as well as p, which 
will be measured through a ratio of vacuum expectation values as described 
below, and ZA(rn — 0), which will be taken from previous measurements by the 
ALPHA group (see also below). 

An equation for the axial vertex. Although not necessary for our direct 
purpose, we write for completeness the equation for the axial vertex parallel to 
the one relating the vector vertex to the quark propagator. Let us recall that 
one could use this relation to express Zy in terms of the bare propagator and 
vertex according to : 

Z v = Z^p 2 )/g^{p 2 ) (17) 

where gy\p 2 ) is the coefficient of the 7^ term in the Lorentz decomposition of 
the vertex. In practice, this relation is not the best suited to measure accurately 
Zy because of many artefacts. 

By taking the derivative of the axial W-T identity with respect to q^, it is 
indeed similarly possible in principle to determine the axial vertex at q — from 
the propagator and the pseudoscalar vertex. 

In the axial case, the expression is more complicated that in the vector 
one; this is due to the last term in eq. ([TTj) which originates in the pseudoscalar 
density, and reflects a non conservation of the axial current. What should be 
stressed is that once more the effect of this term does not vanish in the chiral 
limit if the symmetry is spontaneously broken. 

As in the vector current case, the W-T identity will give the constraints on 
the axial vertex by taking the derivative of eq. (|TTj) with respect to q at q = 0. 
We get : 

ZaTa, = -i-^S- 1 (p) l5 + 2iZ A p-^T 5 (p) (18) 
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This relation again shows that Za is independent of the renormalisation scheme. 
Of course, once more, this will hold up to terms vanishing as inverse powers of 
the cutoff at infinite cutoff, which are called artefacts in the lattice language. It 
must be recalled that, on the lattice, the Ward identity is not exact, but holds 
only up to artefacts, because we work at finite cutoff, and the deviation will be 
found very large in some cases. Although theoretically possible, in practice, it is 
not easy to determine accurately Za from it. This is why we recourse to other 
determinations from the ALPHA group. 



2.3 MOM renormalisation constants 

Useful and very usual specific renormalisation conditions for the Green func- 
tions on the lattice are the MOM ones, considered at some normalisation mo- 
mentum p 2 — fi 2 , originally due to Georgi and Politzer. Although they are not 
really needed for our purpose, it remains useful to explain the connection with 
what precedes, since results are often discussed in terms of the renormalisation 
constants corresponding to this scheme. We lay particular emphasis on the need 
to account for Ward-Takahashi (W-T) identities in handling renormalisation. 
We start from the propagator and set : 

Sr 1 ^ 1 ) = 5 a,b {i1> + TOr) | p 2 =/j 2, (19) 

which means Z^ /OM = Z^(/i 2 ) -1 according to eq. (fT]). Also, it means that 
the renormalised mass is then m^ OM = m([j, 2 ), i.e. it is the mass function at 
p 2 = /i 2 . This is the scheme of Georgi and Politzer. 

For the vertices, we could think of choosing also the standard MOM ones, 
i.e. tree level expressions for p 2 — /i 2 , namely : 

(^W) = i 

(gPW) = i (20) 

This would be satisfactory at zero mass in perturbation. However, we must 
recall that one is not free of choosing the renormalisation of vertices once the 
scheme has been chosen for the propagator. Indeed, the renormalised theory 
must obey the symmetries, and this fact translates itself into W-T identities 
strongly constraining the vertices. 

In the case of the vector current, the W-T constraint still implies that : 



Zy = ^l=Z^M {fi% (21) 
9v (M ) 

where Zy OM is defined by the standard MOM renormalisation condition : 

(^W) = I- (22) 

But this does not work for the axial current. In effect, it is necessary to 
deduce the renormalisation of the axial vertex from the W-T identities in a non 
perturbative treatment, with spontaneous symmetry breaking. From the W-T 
identities, one deduces that Z A ± Zf OM , where Zf OM would be defined, in 
parallel with Z$ IOM , as Z.^ 2 ) / g A 1] (p 2 ) ; in fact Z^(fj?) / ' g [ A 1] (> 2 ) is not even 
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independent of p 2 ; one can expect only that it reaches Za at large p ; then, 
it would be perhaps better to discard this MOM definition, since it could be 
misleading. 

On the reverse, Zp can be defined, consistently with the Ward identities, 
from the standard renormalisation condition, analogously to Zy . Indeed, setting 
m MOM _ tti^ 2 ) and making p 2 = y? in eq. (fl"2")l . one gets : 

W) = 4^ = ^ /OA V), (23) 

showing that Zp(/j, 2 ) as deduced from the W-T identity and MOM conditions 
for the propagator is equal to the one standardly defined directly by the tree 
condition : 

(fl^W) = I- (24) 



2.4 OPE in the chiral limit and the quark condensate 

The quark propagator, like any Green function, can be described by OPE at 
large momenta. In the present case, one deals with a non gauge invariant Green 
function, which implies the potential presence of non gauge invariant local op- 
erators in the OPE. The coefficient of all the operators giving the leading power 
corrections have been calculated by Lavelle et al. [T2]. A great simplification is 
obtained if one takes the chiral limit. Then, for the scalar part of the propa- 
gator or of its inverse, or the quark mass function, the perturbative contribution 
vanishes, and the dominating contribution in the OPE is the quark condensate 
one ; this is in agreement with the fact that both the quark mass function and 
the condensate vanish with restoration of chiral symmetry . We stress this ex- 
ceptional situation where the OPE begins by a power correction. This could 
lead in principle to a complementary determination of the condensate, but as we 
shall see in section [SJ the attempt fails. The relevant quantitative formulae are 
given in this latter section. We just recall the tree level formula proposed a long 
time ago by Politzer (13j , which gives the general structure of the contribution : 

m{p 2 ) - -^r" s (p)4 < W > ( 25 ) 

2.5 A calculation of the condensate through the pseu- 
doscalar correlator 

A very interesting identity has been considered some years ago in several 
works on the lattice, which leads to a possible calculation of the condensate, ad- 
vantageous for our discussion below, subsection l4.2l This identity holds exactly 
when chiral symmetry is preserved by the lattice regularisation, and it has been 
indeed written for overlap fermions , see references |14[|15j : 

< >= - lim m d 4 x< P 5 {x)P 5 (0) > . (26) 

m— >0 J 

as adapted to our specific case Nf = 0. The quantities are defined with the 
"rotated" quark fields, often introduced in overlap calculations to improve Green 
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functions (see the above references). Note that for m ^ both sides retain 
additive divergences proportional to powers of the mass, although the strongest 
one ,1/a 3 , is cancelled by the quark field rotation. On the l.h.s, these divergences 
are due to the limit x — > taken in the propagator. On the r.h.s., the integration 
may imply divergences from coinciding arguments of the two P5 composite fields, 
in addition to the multiplicative ones coming from each P5 field, so that these 
additive divergences cancel between the two sides. In the chiral limit, we have 
not such divergences, since they are cancelled by powers of m. In renormalised 
form, using for the overlap action Zp = Zg = Z.^ 1 , we get : 

< W > R = - lim (m R fd 4 x < (P 5 )r(x)(P 5 ) r (0) >) (27) 

Now, in the renormalised form, it should hold independently of any par- 
ticular choice of action, therefore also with the clover action, up to artefacts. 
Taking duly into account the renormalisation factors, we end, with the bare 
< P 5 (x)P 5 (0) > on the r.h.s. , with : 

< M >R= ~ lim (Z s {Zp/Z s ) 2 m [ d A x < P 5 {x)P 5 (0) >) (28) 
m-R—^O J 

The ratio Zp/Z$ is independent of the renormalisation scheme. Only Z$ de- 
pends on it. The bare < Ps(a;)P5(0) > on the r.h.s. can now be taken from any 
action, including the clover one. 

This identity can then be used on the lattice to calculate the condensate 
rather directly from the pseudoscalar correlator, for instance in the MS scheme. 
Indeed, the r.h.s. can be calculated with only standard logarithmic multiplica- 
tive renormalisations. We need not extract the pion residue, although this can 
be done as in the quoted papers, ending on the GMOR relation. It is useful for 
our purpose to avoid the recourse to the GMOR relation, because we want to 
calculate the condensate in the absence of the Goldstone state. In fact, it can 
be noticed that this direct method is the one proposed a very long time ago by 
Bochicchio et al., the Rome group, [7J. It can be established starting from the 
standard axial W-T identity : 

d M < A^{x)(P 5 ) R >= 2m R < {P 5 ) K {x)(P 5 ) R > +6(x) < ^ > R (29) 

and integrating over x. Note that we duly postulate the Ward identity in renor- 
malised form. It is valid in bare form only up to finite renormalisation factors 
Zp/Zs- This Ward identity has been first established and exploited by D.J. 
Broadhurst |16j . 

This method is advantageous also with respect to the one of calculating 
directly < tpip > through the propagator : it seems to circumvent the problem 
of extracting the power divergences. 



3 Lattice calculations with the clover action 

3.1 The problem of the Wilson term circumvented by the 
study of the vertex 

Of course, versions of the lattice Dirac action have been devised to improve 
the chiral behavior of the Wilson action, like the overlap, domain wall or twisted 
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fermions. The advantage of the Sheikholeslami-Wohlert(SW) or clover fermions 
is that they are relatively easily handled for not too small masses, in contrast 
to these more sophisticated versions. 

To calculate m(p 2 ), the simplest way would seem to extract it directly from 
the lattice propagator, by extracting the scalar part. But this is not practicable 
for the Wilson or clover action, due to the large magnitude of the Wilson term, 
which affects the scalar part. Indeed, it is of order 0(ap 2 ), where a is the lattice 
unit ; this term is purely an artefact, but it cannot be avoided ; not only it 
is large, but moreover it increases like p 2 , while the real, continuum m(p 2 ) is 
decreasing like 1/p 2 in the chiral limit, and otherwise logarithmically. Let us 
remind that the clover action improves the on shell quantities from order 0(a) 
down to 0(a 2 ) [T7], but not the Green functions, and presents the same large 
Wilson term artefact. 

Fortunately, the problem is circumvented by the study of the vertex [3], 
because, in this case, no such an embarrassing artefact is present. Indeed, as 
we shall see, the data on the pseudoscalar vertex, with only a proper treatment 
of hypercubic anisotropy, present the expected, roughly power-like, decreasing 
behavior, see Fig. [TJfor the related m(p 2 ) oc Zp(p 2 ). Moreover, the smallness 
of artefacts is guaranteed by the good superposition of the data at 6.0 and 6.4 
with proper renormalisation. This weakness of artefacts is only an empirical 
fact devoid of explanation. One may suspect that it is connected with the 
"amputation" of external lines, i.e. for instance the tree approximation is exactly 
75 ; nevertheless, this is not a sufficient reason, for the analogous vector vertex 
is found to be still spoiled by large artefacts [MTU], 

The connection between m{p 2 ) and the pseudoscalar vertex is given in eq. (|T5|) 
of the previous section, up to artefacts. Thereafter, for simplicity of notation, we 
denote by m(p 2 ) the combination given by the r.h.s. of this equation, which 
should be equal to the scalar part of the propagator in the continuum limit, but 
which is quite different on usual lattices. We introduce lattice units : 

m(p 2 ) = a- 1 Z A p gf ] (p 2 ) / Z^{p 2 ). (30) 

p is now the dimensionless bare axial quark mass : 

ad v f A = 2pP 5 (31) 

with j A and P5 the bare lattice local axial current and pseudoscalar density 
and a is the lattice spacing. In practice, p is determined in the standard way 
through the ratio of v.e.v.'s : 

' ' E s < P 5 {x,0)P 5 {0) > K > 

Hereafter, we will define k c as the value of k at which p vanishes. All the other 
factors are dimensionless ab initio. Z A is the p independent renormalisation of 
the local axial current. It can be determined by certain W-T identities among 
current correlators j!8) or other methods with specified renormalisation condi- 
tions. The lattice definition of Z^p 2 ) is, as previously : 

MP) = J^Tr [j^S-\p)] /{pf (33) 
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where p M = i sin (ap^). Equation (f3"0"f is the well known formula which has been 
used classically to determine the renormalised quark masses on the lattice [TO] 
at short distance. But here, it is used in the non perturbative regime of 
spontaneous chiral symmetry breaking, and in the chiral limit where m q = . 
In order to avoid an increase of errors on the chiral extrapolation, and since we 
are interested only in the chiral limit of m(p 2 ), we multiply by Za only after 
having taken the chiral limit of the remaining factors. Now, the Za§ in the 
chiral limit are accurately known by the work of the Alpha group |20j . and we 
borrow their central values Za(k c ) . 

3.2 The treatment of the raw lattice data 

The basic data, i.e. the quark propagator in the various configurations, are 
the same as already used in [9UTD] to study Z^{p 2 ). We have at hand simulations 
at Nf = with the Wilson gauge action and the SW clover action in Landau 
(i.e. Lorentz) gauge, on a series of lattices given in order of decreasing physical 
volume : 6.0, 24 4 L = 12.2 GeV' 1 ; 6.0, 16 4 , L = 8.14 GeV 1 ; 6.4, 24 4 , L = 
6.56 GeV' 1 ; 6.6, 24 4 , L = 5.06 GeV' 1 ; 6.4, 16 4 , L = 4.37 GeV' 1 ; 6.8, 24 4 , 
L = 3.93 GeV . We have considered the inversion of the Dirac operator at 
five kappa values in each case ; namely at /3 = 6.0, we choose re ranging from 
0.1310 to 0.1346, corresponding to a large range of quark masses, so as to allow 
a reasonable chiral extrapolation ; the values at the other /3's are chosen to 
correspond approximately to the same bare masses in physical units, i.e. in 
terms of m q = ^(l/re-l/re c ), m q = 0.233,0.154,0.104,0.054,0.0324 GeV. 

We first calculate the product pg^\p 2 ) / Z^{p 2 ) at each re, for a given (3 
and volume 16 4 , 24 4 , and we correct this quantity for the hypercubic artefacts 
according to the same method used for the propagator vector part Z^,(p 2 ) (see 
discussion below). Then we take the chiral limit re —> re c according to eq. (|30l) . to 
obtain the chiral limit of the quark mass function. The chiral limit is obtained by 
a fit in function of re. The factor Za(k c ), which has a trivial effect is introduced 
as an additional fixed factor, not affecting the p dependence and the essential 
conclusions. A three parameters fit in re is possible with five re values, and it 
gives a sizeably better fit than with two parameters, pointing to a significant 
curvature in m q . The coefficient of the 0(m 2 ) term is strongly negative. 

3.3 Discussion of artefacts 

We have given a very detailed discussion of artefacts in our previous papers 
on the quark propagator |9ll0j . We address the reader to these papers. However, 
it happens that some aspects of the discussion are crucial here, so that we give 
a new discussion for the relevant points. 

3.3.1 Discretisation artefacts 

As to discretisation artefacts, let us recall that we can classify them into 
two categories according to their behaviour under 0(4) : either they are not 
invariant under 0(4), which corresponds to the hypercubic artefacts, and we 
may determine them by using the various orbits ; or they are invariant and 
they can be extracted by using various /3's. It happens that in the present case, 
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contrarily to the vector part of the propagator, both hypercubic artefacts and 
0(4) invariant artefacts seem small. 

Note that we have treated the hypercubic artefacts by our systematic method 
of extrapolation, explained in several places, see for instance [2T] 3 . With this 
method, we obtain data where almost any anisotropy has been eliminated, see 
Fig. Q]for an example, with the chiral extrapolation performed. 
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Figure 1: The chiral extrapolations at (3 — 6.0, with 24 4 and 16 4 lattices, after 
elimination of hypercubic artefacts. A very regular behaviour is observed. The spacing 
between the two curves increases in the IR, which can be interpreted as a finite size 
effect. 



With such smoothed data, it is possible to perform very good analytical fits 
based on theoretical considerations in the continuum, with very low \ 2 ■ Possibly, 
the fit should also include terms accounting for 0(4) symmetric artefacts, to be 
determined by considering several j3's. In the present case, we even do not 
require important terms of this sort. We have a good superposition of the chiral 
extrapolations at 6.0 and 6.4 with the 24 4 lattice, on their common range of 
momenta, within statistical errors, see Fig[2] However, a further, more detailed 
study, at large p, with multiplication by p 2 , reveals the possibility of a small 
0(4) symmetric artefact (see the end of the discussion of OPE, section [3] and 
Fig. E}. 



On the other hand, at 6.6 and 6.8 for the same lattice size, except for the 

3 This technique was initially devised by C. Roiesnel 
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Figure 2: The chiral extrapolations at /3 — 6.0 and /3 = 6.4 at 24 superpose very well 
on their common range of physical momenta. There remains however a small finite 
size effect. 



first few points (where we have most probably a finite size artefact, see below), 
the chiral extrapolations superpose very well, but at a zero value, therefore 
completely different from the previous case, see Fig. [3] Obviously, this is not a 
discretisation problem. We show in the next section that there is a discontinuity 
in function of the physical volume. 

Having thus explained our own results, we give some comments on a recipe 
for treating discretisation artefacts , which leads to surprisingly different results 
in the case of the pseudoscalar vertex [52], in the chiral limit. In addition 
to an usual democratic selection, reference (22) reads the continuum p M as cor- 
responding, on the lattice, not to the lattice p^, but rather to a trigonometric 
expression, sin(ap^) / 'a, differing from p^ by 0(a 2 ) terms. As far as the residue 
of the pseudoscalar vertex is concerned, our systematic method for eliminating 
hypercubic artefacts happens to give a result close to the democratic selection 
if one reads p 2 in the democratic method as the S/ij? 2 of the lattice. However a 
large discrepancy appears when [22] identifies p 2 to S /J (sm(ap M ) /a) 2 , for the 
very asymmetric lattices used there (for example, 16 3 x 52 or 24 3 , 64). Why the 
effect can be large can be easily understood : the residue behaves roughly as 1/p 2 
in physical terms, and then, with the sine recipe, the curve appears much lower 
at large p, because the ratio of the two curves is roughly E(sm(ap p )/a) 2 /E AI p 2 I ; 
one finds differences as large as 50%, for some asymmetric lattices at the largest 
momenta. Then, at large p, the recipe gives a notably smaller result than we 
find. Now, the question is : what is the correct answer ? 

Our method gives a definite answer, by exploiting various orbits. As to the 
recipe, it could seem to be somewhat justified in some particular places, for 
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Figure 3: The chiral extrapolations at ft = 6.0 and f3 — 6.4 superpose very well, with 
a relatively large value, corresponding to spontaneous chiral symmetry breaking. On 
the other hand, the chiral extrapolations at ft = 6.6 and /3 = 6.8 also superpose very 
well, but with a zero value, corresponding to unbroken chiral symmetry. The huge 
difference between the two sets is a physical volume effect, as explained in section [4] 



free quarks, by comparing the continuum and lattice explicit expressions of the 
Green function : for example, the vector part of the inverse quark propagator, 
which is found to be 1X7^ sin^ap^/a on the lattice, instead of «S7 M p^ in the 
continuum. Nevertheless, this has been taken into account by the standard 
lattice definition of Z^(p 2 ), eq. (|33l) . which gives exactly one for free quark, like 
in the continuum ; now, when one considers not the free Green functions, but 
the non trivial p 2 dependence of the invariants like Z^ip 2 ) or g^\p 2 ) due to 
the interaction, the recipe has no theoretical justification, and one can even 
doubt that there is a universal empirical recipe to reduce UV artefacts 4 . 



3.3.2 Lattice finite size artefacts 

The last type of artefacts is what will be usually termed as finite size arte- 
facts, affecting the first few points in momentum space, i.e. the one with lowest 
momentum number. Let us stress that it is quite distinct from the effect of 
the physical volume, which we study in the next section, and which extends 
all over the range of available momenta and shows a discontinuity. 

As to these finite size artefacts, we find that they are present. In fact, we 
observe that whatever the number of sites and the lattice unit, the first five 
points are always enhanced. We show this by considering first the cases where 
m(p 2 ) is found to be very close to zero, Fig. @]: although the /3's are different, 

4 The Adelaide group has indeed observed, in the case of the overlap action, and studying 
the scaling behaviour, that two different such recipes should be used for the vector and scalar 
part of the propagator 
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we can see that the value is extremely small (10 -4 GeV !!!) for 6.6 and 6.8 
except for the first five points ; in this case, they differ from each other, which 
shows that it is an artefact. The data at 16 4 , which are larger, although still 
very small (10~ 3 GeV !!!) except at the three first points, suggest the same 
interpretation, but also suggest that the finite size effects are notably larger at 
16 4 . 
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Figure 4: The cases /? = 6.6 or /3 = 6.8 at 24 4 , and /3 = 6.4 at 16 4 , corresponding to 
the restoration of chiral symmetry. One sees IR finite size effects, especially large in 
the latter case. 



Then we observe that there is a similar effect in the other case m(p 2 ) ^ 0, 
by comparing the 24 4 and the 16 4 lattice at the same 6.0, Fig. [TJ or (3 = 6.0 
and P = 6.4 with size 24 4 , Fig. H 

4 Discontinuity of chiral extrapolations in func- 
tion of the physical volume 

4.1 Discontinuity of the quark mass function 

We observe two very distinct cases, already seen in Fig. [3]for the 24 4 lattices : 
a) at (3 = 6.0 and 24 4 or 16 4 , as well as at f3 = 6.4 and 24 4 , we obtain a non 
zero chiral limit, which superposes very well where possible, in the three cases, 
with large values at small p, i.e. of order of 200 MeV at p ~ 0.5 GeV. This 
typically exhibits the behavior of a spontaneous breaking of chiral symmetry. 

The presence of a curvature in function of m q seems to explain why our 
chiral limit is smaller than in a previous calculation [3], where only three k's 
were available : for illustration we were getting around m(2 GeV) ~ 0.034 GeV, 
while we obtain now m(2 GeV) ~ 0.022 GeV. 
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b) at /3 = 6.6 and 24 4 , and at /3 = 6.8 and 24 4 , as well as at f3 = 6.4 
and 16 4 , we obtain on the contrary much smaller values on the whole range 
of momenta ; moreover, if we recalculate k c by defining it through p{n c ) = 
in each corresponding finite volume, it is a bit different from the standard one 
determined in very large volumes, and then the chiral limit is very small and 
physically not significant 5 , except for the low momenta p < 3 GeV of the 
6.4, 16 4 lattice, where there is a value significantly different from zero, although 
small, and monotonously decreasing with p. For these latter points, the most 
natural explanation seems a volume artefact on the four or five points, which is 
also present in the other cases, with a different magnitude, as we have argued. 
On the whole, the case b) seems typically a situation of restoration of chiral 
symmetry. 

Now, the remarkable fact is that chiral symmetry "seems" to be "restored" 
rather abruptly for volumes smaller than a certain physical value L c ~ 6 GeV^ 1 . 
in the sense that standard chiral extrapolation by a low polynomial in m q yields 
m chl rai{p 2 ) above L = L c and m chlral (p 2 ) = below, where m chlra i denotes 
the chiral limit. Indeed, the relevant parameter distinguishing case a) and case 
b) seems to be the physical volume. The physical length of the lattice is respec- 
tively L = 12.2, 8.13, 6.55 GeV^ 1 for the case a), and L = 5.05, 4.37, 3.93 GeV' 1 
for the case b), the order being the same as above. The separating length is 
then around L c ~ 5.8 GeV -1 = 1.2 fm (we choose the middle between the 
two lengths). Moreover, we speak of a discontinuity and not simply of a tran- 
sition : it is because the last volume which presents "symmetry breaking" is 
certainly larger than the first one which presents "symmetry restoration" , but 
the two volumes are not very different : (3 — 6.4, 24 4 ,L = 6.56 GeV^ 1 against 
f3 = 6.6,24 4 ,i = 5.06 GeV~ . Moreover, on each side of the discontinuity, the 
results for the order parameter m c hirai(p 2 ) are very similar for the three lattices 
of case a), and very similar for the three lattices of case b). In the case of 
spontaneous breaking, this is especially striking, because the volumes extend 
over a large range ; the largest volume, f3 = 6.0, 24 4 can be considered as a 
relatively large volume, V = 2.2 10 4 GeV~ 4 but the smallest one /? = 6.4, 24 4 
which is ten times smaller, V = 1.85 10 3 GeV"* ; and the latter volume is much 
closer to the /3 — 6.6, 24 4 volume V = 6.55 10 2 GeV~ 4 for which symmetry is 
manifestly restored. This means a discontinuity around V = 1.25 10 3 GeV~ A 
(middle point) or, in length, around 5.8 GeV -1 . 

In [TU], we studied the chiral symmetry breaking through the overlap r a/IV 
as a possible indicator of spontaneous symmetry breaking : it should differ from 
1 due a Goldstone contribution at small p ; in fact, we saw an effect at small 
p presumably coming from the Goldstone, at least at j3 — 6.0. However, at 
small volumes corresponding to 6.6, 6.8, the available momenta were too large 
to test whether this Goldstone effect disappears or not. Now, the study of the 
pseudoscalar vertex (or quark mass function), for which the Goldstone effect is 
much larger when present, gives a clear answer : it shows a striking vanishing 
of the Goldstone effect at 6.6, 6.8. 

A typical feature is that this transition affects m c hirai(p 2 ) over the whole 
momentum range simultaneously, i.e. when passing the same critical length. 
Of course, there are also, as we have shown in subsubsection 13.3.21 finite size 

5 In view of this observation, we have recalculated also the case a) with the same prescription 
for k c . 
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artefacts which somewhat enhance the smallest momenta. But they superpose 
on top of a very clear discontinuity with respect to the physical volume. And, 
even where enhanced by this volume artefact, i.e. at small momenta, the curves 
of case b) lie much lower than those of case a) : admittedly, the 6.4, 16 4 is also 
lower, but less than the others at the very first few points. 

To reinforce our conviction that we are indeed facing a notable phenomenon, 
we propose in the next subsection two other similar and striking observations, 
which display the same critical length. 

4.2 Discontinuity of the pion mass and of the condensate 
value 

We extract these two observations from the behaviour of the pseudoscalar 
correlator, which we have calculated together with the vertex on the same lat- 
tices with the same four f3 values. From it, we can calculate, 

i) the pion mass, according to the well known method using the< P 5 (x)P 5 (Q) > 
correlator integrated over space, at large t, 

but also 

ii) the quark condensate, through the method explained above (subsection [53]), 
i.e. through the chiral extrapolation of — m J d A x < P^(x)P5(0) >. 

For both quantities, on performing a low order polynomial extrapolation in 
m q = Y/2a(^-/ K — V K c)j we observe exactly the same type of discontinuity as 
for m(p 2 ), i.e. everything goes as if one had a phase transition at small physical 
volume. The results are displayed in the table below and graphically in figure 
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p 


Size (Length) 
( GeV" 1 ) 


(GeV 2 ) 


(GeV 3 ) 




6.0 


24 4 (12.2) 


(1.54 ± 1.58)10^ 


(-44.3±7.4)10- 3 


a) 

"broken 
symmetry" 


16 4 (8.14) 


(4.05±2.59)10~ 2 ~ 


(-50.7±7.9)10- 3 


6.4 


24 4 (6.56) 


(-5.1 ± 4.6)10-* 


(-43.5±7.8)10" d 


6.6 


24 4 (5.06) 


(0.66 ±0.115) 


(-7.29 ± 14.3) 10" 3 


b) 

"restored 
symmetry" 


6.4 


16 4 (4.37) 


(0.615 ±0.21) 


(-24.5 ± 12.1)10 -3 * 


6.8 


24 4 (3.93) 


(1.3 ±0.22) 


(-8.65 ±28.8)10^ 



We draw the following conclusion from this table : 

i) the chiral extrapolation gives a much smaller m 2 in case a) than in case 
b), by two orders of magnitude. This striking difference confirms the advocated 
discontinuity. In fact, it is almost compatible with 0, m 2 ~ 0, in case a) - 
and this is in agreement with spontaneous symmetry breaking. For the largest 
volume, it is fully compatible with 0. On the other hand, there is no Goldstone 
boson in case b), but a heavy meson. A complementary criterion for the study 
of chiral restoration would be the presence of a scalar meson degenerate with 
the pseudoscalar one. 

ii) the chiral extrapolation of m J cftx < Ps(x)P5(0) > gives something close 
to the expected value for the quark condensate, or larger than it for the case a). 

The central values of the corresponding renormalized condensates in MS at 
2 GeV are estimated to be, either by using values of z^*~ non perturba- 
tively measured or those obtained through Zp ~ , and using also Zp/Z$ 
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from Ward identities :< inp > TT g (2 GeV) = -(2.2±0.4) 10" 2 GeV 3 at 6.0, 24 4 , 

< W >ls ( 2 GeV ) = ~( 2 - 5 ± °- 4 ) 10_2 Gey3 at 6 - ' lfi4 > m a q uite encour- 
aging agreement with standard values from the GMOR relation at Nj = 0. 
Indeed, from [35] : 

< tpi/j >jjg (2 GeV) = -(273 ± 19 MeVf = (2.0 ± 0.5) 10~ 2 GeV 3 (34) 

On the other hand, the renormalised value < tjji/j >jfg (2 GeV) — —(3.2 ± 
0.6) 10~ 2 GeV 3 obtained at 6.4, 24 4 is somewhat too large, but not far from 
compatibility with eq. 

On the contrary, for the case b), it is smaller, and compatible with zero for 
two lattice out of three. The exception is the 6.4, 16 , which gives an unexpected 
intermediate value. So the conclusion would be equally striking , if not for the 
exception of this small volume lattice 6.4, 16 4 . 

This exception does not correspond to an intermediate physical volume (this 
would invalidate our advocated conclusion of a discontinuity as function of the 
physical volume). Rather, it is probably due to a finite size artefact, connected 
with the similar observation for m(p 2 ). At the same 6.4, 16 4 , m(p 2 ) shows 
rather high points at the first momenta, much higher than those for the others 
in case b). Such finite size artefacts seem also to affect the case a) ; the renor- 
malised value of the condensate at 6.0, 16 4 is somewhat larger that the one at 
the largest volume 6.0, 24 4 ; simultaneously, the pion mass is somewhat larger. 
It is difficult to give a complete rational explanation of the size of such artefacts. 
Nevertheless, it is doubtless that finite size artefacts are superimposed on the 
basic physical volume effect, and can be distinguished from it because they do 
not follow the same rationale. 

4.3 Comments on the unexpected discontinuity of chiral 
extrapolations in function of the physical volume 

4.3.1 A phase transition ? 

At this point, it is important to enter into some warnings, to avoid possible 
confusions. 

i) What is as expected. We are well aware that chiral symmetry is not ex- 
pected to undergo spontaneous breakdown in finite volumes : order parameters 
are expected to vanish anyway when m q — > 0. But, if the volume is suf- 
ficiently large, this may happen at only very small quark masses, which we 
cannot consider with Wilson-type actions. Then, at moderate masses reached 
by Wilson-type actions, things, as is well known, may look as in the infinite 
volume case : order parameters like the mass function seem to tend to the same 
non zero value as in the infinite volume, if one performs a naive chiral extrap- 
olation with a low order polynomial in m q . Of course, m(p 2 ) as function of 
Tfiq should present a bending downwards with respect to this extrapolation, if 
one were able to go to smaller masses. To know whether the volume is suf- 
ficiently large, a necessary condition is to check whether the "chiral" limit is 
stable against variation of the volume, which we verify. Let us call this first 
situation "infinite- volume-like" . It seems to correspond to our case a). 

On the other hand, if the volume is sufficiently small, it is natural to expect 
that, simply extrapolating from the moderate masses, one gets already a van- 
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ishing of order parameters. Let us denote this second situation as the one of 
"chiral symmetry restoration". It seems to correspond to our case b). 

ii) What is surprising. This dichotomic presentation is however an oversim- 
plification, according to the common ideas. One would expect a continuous 
transition between the two situations when one decreases the volume ; the 
chiral extrapolation would be expected to deviate progressively from the truly 
infinite volume value, and decrease downwards to zero. In the same vein, the 
Goldstone boson would be expected to acquire progressively a mass. Let us 
recall for instance the finding in the so called e regime for the condensate : the 
ratio to the infinite volume limit deviates from 1 by a function of z — m V, 
therefore, at fixed mass, it is a continuous function of the volume 

In face of such expectations, our lattice analysi show on the contrary an 
abrupt discontinuity between the two cases a) and b) in a rather narrow win- 
dow of physical lattice length or volume. Indeed, our low order polynomial ex- 
trapolations exhibits this striking discontinuity for m c hirai{p 2 ), which drops 
suddenly to zero. And the Goldstone boson acquires abruptly a very large mass. 
This discontinuity would suggest speaking of a phase transition. In fact, this 
would be a too strong statement, since principles seem to be against such a 
conclusion for finite volumes, and also, since we are not really experimenting 
the chiral limit, with truly very small masses, but instead performing only an 
extrapolation. Rather, one should speak of a sharp transition to a new regime 
of chiral extrapolation : 

-above the critical length, the naive chiral extrapolation picks the quantities 
corresponding to infinite volume and spontaneous breaking, 

-while thereunder it picks the ones corresponding to chiral symmetry restora- 
tion. 

This behaviour is illustrated in figure [5] 
4.3.2 Origin of the discontinuity 

Our observation deserves obviously understanding. We are tempted to assume a 
connection with the observation of Neuberger and Narayanan [5] of a restoration 
of chiral symmetry below some critical length. Admittedly, they work in the 
Nc — > oo limit and they then expect a true phase transition at finite volume. 
However, their critical length is close to ours : they find L cr i t around 1 fm — 
5 GeV -1 . 

One can also think of a connection with the "finite temperature transition" , 
which affects both the confining properties and the chiral symmetry and is ob- 
served at large spatial volume and finite time interval 1/T. Indeed the transition 
occurs around T c — 0.270 GeV, therefore 1/T C ~ 3.7 Gey 1 , mot too far from 
our L c = 5.8 GeV^ 1 . However, our data do not correspond to this situation of 
highly asymmetric lattices which would allow to study such questions. 

Let us emphasize that the transition we observe, concerning chiral symmetry, 
is not a universal fact. Not all condensates vanish at small volumes. In our 
analyses of the gluon, ghost and vector part of the quark propagator, the non 
perturbative v.e.v < A 2 > was found the same, consistently, at all the four 
volumes j3 = 6.0, 6.4, 6.6, 6.8, 24 4 , the same volumes as we use now. < A 2 >, 
which does not seem to be an order parameter for some symmetry breaking, 
does not collapse in small volumes. Note that we are not yet in a situation 
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Figure 5: m% (left) and the quark condensate (right) as functions of the lattice size. 
The leftmost part of each graph shows the restoration of chiral symmetry at small 
volume while the symmetry is broken in the rightmost one. The horizontal lines are 
indicative of the average value of the measured quantity in each region. Note that for 
in the symmetry breaking phase the errors are very small and can hardly be seen 
on the figure. 



of such a small volume that everything should be perturbative simply because 
power corrections would be negligible ( "femto- universe" ) . 

5 The OPE 

5.1 Failure of OPE with the Wilson coefficients at low 
order 

Let us now concentrate on the case a), i.e. the one where the chiral limit mass 
function is large. From now on, we shall consider only renormaliscd quantities, 
in the MS scheme, and skip any subindex meant to recall the renormalisation : 
we quote only the subtraction point. Let us first recall the puzzle underlined 
in [3], where one studies the results of lattice data by the QCDSF group at 
j3 = 6.0. If we consider the OPE with tree level coefficient for the condensate, 
we have in the chiral limit the Politzer formula quoted above [13] : 

m(p 2 ) = -^-a s {p)\ <Tpip> (35) 
3 p A 

Note the remarkable fact that the expansion begins with the power correction. 
The purely perturbative contribution vanishes in the chiral limit, since it is 
proportional m q . 

We observe immediately an enormous discrepancy between both sides 
at the momenta usually considered ; indeed, with the estimate of < >= 
-[(0.267±5±15) GeV] 3 = -(1.9 ± 0.35)10" 2 GeV 3 from the mass of the pion 
and GMOR formula with overlap fermions ( [IS]), or — (2.±0.5)10~ 2 GeV 3 from 
clover fermions, as given above [22], see eq. ; taking p — 2 GeV, a s (p) ~ 0.3 
(at one loop, with Aj-j-g = 0.240 GeV), we find a r.h.s. around 5 10~ 3 GeV, 
much smaller, by a factor 4, than the value we find from the lattice result for 
the l.h.s. : miatt(2 GeV) ~ 2 10 -2 GeV .Note that there is an habit to present 
condensates through cubic roots, which hides the discrepancies. We want to 
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avoid it. Le us recall that 4 1 / 3 ~ 1.6. Note also that, in the initial paper [3J, 
the discrepancy was found still larger, around 10. The reason is twofold : first 
we adopted as reference the standard QCD sum rule value of the condensate, 

< tfnf) >= -(0.225 GeV) 3 = -0.0114 GeV 3 ,which is twice smaller ; second, the 
lattice value of m(p 2 ) that we estimated was larger, as explained in the beginning 
of section^ because we could do only a linear extrapolation. Nevertheless, the 
discrepancy remains huge. 

Working at higher momentum would not help much in this respect : us- 
ing the (3 = 6.4, 24 4 lattice, we have momenta up to more than 5 GeV, but 
nevertheless the discrepancy is not much smaller : the prediction at tree level, 
with a s (5.7 GeV) ~ 0.18 m(5 GeV) ~ ^0.18^0.0176 GeV ~ 4 10~ 4 GeV, 
against mi at t{5.7 GeV) — (1.15 ± 0.67) 10~ 3 according to our lattice measure- 
ments, therefore there is a factor about three of discrepancy. Note that here we 
choose to calculate a s (p) at one loop for simplicity in this first discussion. To 
be more quantitative, it seems that we should adopt the best possible approx- 
imation for a s , i.e. including ^3. This is what we do in the rest of the section . 
At Nf = 0, the resulting values of a s are sizeably smaller, which reinforces the 
problem. 

We can include, as done already in [3J, the one-loop correction to the co- 
efficient, which has been calculated a long time before by Pascual and de 
Rafael [23] ; we set in their formula a = for the Lorentz gauge, iV c = 3, 
and for the renormalisation point, /1 = p in their notation ; we also take their 
p 2 as minus the Euclidean p 2 of the lattice ; whence : 

«V) = ~ot a (p)\ <H>> (P) fl + 6.1875^M) (36) 
3 jr \ it J 

We can in addition make a fit on the whole range of our data. We take into 
account the evolution of the condensate from p down to the reference point 

2 GeV where we want to determine the condensate, by using the formula for 
the evolution of the MS quark mass, which is just the inverse of the one for 

< 'tpip > (p) : 

< > (p)/ < fa > (2 GeV) = ((a s (2 GeV)^f llx 
(1 + 0.687328 (a s {2 GeV)/n) + 
1.51211 [a s {2 GeV )/tt) 2 + 4.05787 {a s (2 GeV)/-*) 3 ))/ 
((a s (p)/n)^ n (1 + 0.687328 {a s {p)/ir) + 
1.51211 (a s (p)/n) 2 + 4.05787 {a s {p)/nf)) (37) 

to the corresponding order. We have checked this formula by calculating it from 
both the expressions in 25 and (3SJ[55] respectively. It coincides exactly with 
the one given in [S]. 

In the common range of momenta at /3 = 6.0 and j3 = 6.4, 24 4 , p = 2 — 

3 GeV^, cf Fig. [5J we find a good superposition of the lattice curves for m(p 2 ) 
at 6.4 and 6.0, and a very consistent fit where the "condensate" at 2 GeV in 
WS scheme would have the fitted value -(0.062 ± 0.01) GeV 3 at /3 = 6.0, 
-(0.056 ± 0.009) Ge^ 3 at j3 = 6.4, therefore around 3 or 3.5 times the actual 
value. This is still clearly unacceptable. The situation is only slightly better if 
we use larger momenta, available only with the 6.4 lattices ; in this case, using 



22 



the full range of momenta, we get -(0.05 ± 0.01) Ge^ 3 6 . 

Obviously, the discrepancy is so huge that we have to find a non trivial ex- 
planation. We can think at once of two types of explanation of the discrepancy, 
preserving the OPE : -either the above condensate value, taken from the pion 
mass, has something basically wrong ; it is hard to believe ; note also that many 
successful considerations have been based on roughly this value, and one would 
have to reconsider a whole sector of particle physics ; -or it is the perturba- 
tive expansion of the coefficient which is strongly modified by higher orders, or 
even, which is simply not valid, at least at available momenta. In [3], the second 
explanation is suggested, i.e. it was suggested that the perturbative expansion 
may merely break down at available lattice momenta, since already the one-loop 
correction is found as large as 50% of the tree level at 2 GeV. It could mean 
that either there is merely no sort of convergence or that the OPE is practically 
useless. Equivalently, it could mean that we have to reach very high momenta 
for the known, low order, perturbative expansion to be valid. However, since 
the decrease of a s is very slow, we would require prohibitively large momenta 
to test the idea on the lattice. 

5.2 The lesson from very high order calculations 

Fortunately, an impressive progress has been performed in the calculation In- 
corporating the higher order terms recently computed by K. Chetyrkin and A. 
Maier 7 gives (cf eq. (HIT) in the appendix) : 

m(p 2 ) ~ -An/?>a s /p 2 < ipip > (p) 

(1 + 6.1875(a s /7r) + 52.9495(a s /7r) 2 + 564.8284(a s /7r) 3 ) (38) 

We write ~ because there is in addition a correction of the same formal order 
from the quark self-energy, S(p), beginning at two loops, which is found to 
be very small(l%). The positive and steadily increasing coefficients are much 
larger than in the purely perturbative series for the propagator(see for instance 
eq. 1461) . It do much in convincing one that the perturbative expansion of the 
Wilson coefficient is not well behaved. For p = 2GeV^, a s ~ 0.2 one gets : 

m((2 GeV) 2 ) ~ 4tt/3 0.2/(2 GeV) 2 (- < ^ > (2 GeV)) 

(1 + 0.394 + 0.215 + 0.146) (39) 

On one hand, it shows that the series is meaningless at 2 GeV. And it implies 
that claims to recover the value of the condensate from data around 2 3GeV 
with only the low orders cannot be justified. On the other hand, it suggests 
that the discrepancy we have observed at lower orders between the fitted and 
the real value of < ipip > could be explained by the strength of high order 
radiative corrections. 

6 Note that even introducing a strong additional 1/p 4 term and with a different treat- 
ment of discretisation effects, the condensate in |22| remains sizeably larger than required : 
-(0.312) 3 GeV 3 ~ -0.03 GeV 3 ; without the 1/p 4 term, it rises to a still larger value, larger 
by a factor (. 792/0. 721) 3 = 1.33, whence -0.04 GeV 3 , i.e. a factor 2 of discrepancy. 

7 We thank them for communicating their work to us prior to publication 
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Working at higher p would presumably not suffice, even assuming some 
convergence of the series. With a s (5 GeV)/ir ~ 0.147/7T, the behaviour of 
the series is substantially improved : 

m((5 GeV) 2 ) ~ -4tt/3 0.147/(5 GeP) 2 (- < ^ > (5 GeV)) 

(1 + 0.289 + 0.116 + 0.058) (40) 

and, as can be seen on figure ^ the inclusion of the higher terms results in 
reducing the gap between the lattice data and the OPE estimation by a factor 
of 2. But is it is unable to reconcile them. Equivalently, the fit gives a ficti- 
tious value for the condensate of —(4. ± 0.6)10~ 2 GeV 3 , twice too large. To 
summarize, not even on our 6.4, 24 lattice, where momenta run up to more 
than 5 GeV , can we apply safely OPE. In the case of Z % p{p 2 ) [TU], we could 
use the 6.6, 6.8 lattices, which provide momenta up to 10 GeV, to study the 
< A 2 > condensate; this is not possible for the quark condensate in the present 
case, because of the collapsing to zero of the quark mass function at those small 
volumes. 

What is encouraging is that the terms of the series are all positive, and, if 
this behaviour persists, it makes understandable that the real curve lies above 
the OPE contribution presently calculated. 




2 4 p(GeV) 6 

Figure 6: The lattice data xp 2 at 6.0 and 6.4, compared with the OPE expectations at 
tree level and at three loops, with an upper value of the condensate (2.25 10 -2 GeP 3 ) 
estimated by combination of overlap and clover estimates in the GMOR method, 
showing the gap between the lattice data and the OPE expectations. 



5.3 Remaining discretisation artefacts ? 

Of course, one should consider the possibility that there remain discretisation 
artefacts. We feel that in spite of the good superposition of 6.0 and 6.4, there is 
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some positive 0(4) invariant artefact (e.g. +a 2 p 2 ), but presumably much too 
small to explain the discrepancy. Signals of a small artefact (see Figure [5] are 1) 
6.0 and 6.4 differ somewhat around the endpoint of 6.0, p ~ 3 GeV (see Fig. [5]) ; 
2) the central curve of the 6.4 points seems too flat beyond 4 GeV, although our 
errors are too large to ascertain this statement. About point 2) : the theoretical 
perturbative curve is falling more and more rapidly with increasing order, as seen 
in FigurelT] One may then expect the exact continuum curve to fall still more 
rapidly, and not to be flat. The flatness of the lattice data could be explained 
by a positive artefact of the type a p 2 or a 2 p 2 , corresponding to an artefact 
Sm(p 2 ) oc a or a 2 . The magnitude of this artefact is however strongly limited 
by the region around p ~ 3 GeV, where both 6.0, 6.4 points are present and give 
close central values of p 2 m(p 2 ) (difference is around 0.4 10 -2 GeV 3 against 
a total magnitude of p 2 m{p 2 ) 5.4 10" 2 GeV 3 at p ~ 3 GeV,4.2 10~ 2 GeV 3 
at p ~ 5.7 GeV) but of course it would be better to have smaller errors to 
strengthen this conclusion 8 . 
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4 4.5 p(GeV) 5 5.5 

Figure 7: The fits of p m(p ) with increasing number of loops in the Wilson coefficient 
of the condensate. The fits are performed over the largest momenta available (6.4), to 
favor as much as possible the convergence the OPE, showing the increasing, although 
slowly evolving, slope. 



6 Conclusions 

General approach to the quark mass function The method of circumventing the 
Wilson term artefact which affects the clover propagator quark mass function, 

8 In fact, the initial errors on the pseudoscalar vertex (Zp ) are very small at large p, but 
the extraction of the residue increases them much. 



25 



by passing through the pseudo-scalar vertex, seems efficient, and yields a result 
no too far from the continuum. 

OPE . The present study give useful indications for the extraction of OPE 
power corrections, at least for elementary Green functions. In the regime of 
spontaneous chiral symmetry breaking, the quark mass function is in principle 
an exceptionally favorable case for the study of the OPE power corrections, 
since 

1) the 1/p 2 quark condensate contribution is the leading order in the chiral 
limit, 

2) the value of the condensate is well known by other methods. 

In fact, it is the best place for a lattice, nonperturbative, measurement of 
a Wilson coefficient. In practice, the present study strengthens the previous 
conclusion that an accurate OPE analysis is a very difficult task. 

The present case illustrates one of the difficulties, which will also appear 
in other cases according to the calculations of Chetyrkin and Maier : the per- 
turbative series giving the Wilson coefficient of the power term converges at 
best very slowly ; it is difficult to work at sufficiently large momenta to have 
it converge better. In the present work, at the largest momenta, even with the 
huge improvement provided by the recent work of Chetyrkin and Maier (three 
loops), a discrepancy of a factor 2 seems to remain with the actual value of the 
condensate ; we find a fictitious - < tjjip > (2 GeV) ~ (4 ± 0.6) 1(T 2 GeV 3 . 
Equivalently, one can say that, practically, the OPE even with the very high 
orders calculated can work only at very high p. Effects of similar magnitudes 
are found by Chetyrkin and Maier in several other cases, concerning the A 2 
vacuum condensate contribution to propagators. The situation about the OPE 
of composite operators is not known. 

This causes a difficulty for the program of renormalisation through the lat- 
tice measurement of MOM renormalisation constants. This is especially true 
if we require a high precision of order of the "percent", as often considered. 
Indeed, let us recall that the final aim of the method is to extract from the mea- 
sured renormalisation constants their purely perturbative part. Of course, 
in principle, we could work at a very large p to kill completely the power cor- 
rections, and at a sufficiently large cutoff, in order that this large p should not 
be affected by appreciable UV artefacts. However, from our findings, power 
corrections due to < A 2 > can be guaranteed to be 1% only around p ~ 10 
GeV. Then the task would be obviously very difficult. 

The alternative method is to perform OPE fits at available momenta. Such 
fits must include first the known purely perturbative part, multiplied by the un- 
known Z, second the power corrections, which are predicted to be important at 
usual momenta, but also, as a third contribution, discretisation terms invariant 
under 0(4), which should be important as well at this degree of accuracy, and 
whose magnitude is unknown. The method is then to exploit the different be- 
haviour of the various types of terms, as function of both p and /?. This method 
of extraction of Z is then also difficult, when the form of the discretisation arte- 
facts (as we have shown for the vector part of the propagator) as well as the 
functional form of the Wilson coefficient, which depends on high orders, are not 
known accurately. Moreover, we have not an independent knowledge of the A 2 
condensate value in contrast to the quark condensate. Let us stress that the 
difficulty of OPE is presently manifest only in elementary, non gauge invariant, 
Green functions, but these are precisely the ones used for MOM normalisation 
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conditions. 

"Restoration of chiral symmetry" . On the other hand, the study has revealed 
the unexpected and remarkable feature, that there is a general discontinuity of 
chiral properties, as obtained by the standard chiral extrapolation methods : 
when the physical volume passes through some critical value, chiral symmetry 
seems to be restored abruptly although we do not expect a real phase transition. 
The discontinuity affects simultaneously the quark mass function, the pion mass 
and the chiral condensate. Since the critical volume coincides with the one 
found for the phase transition to symmetry restoration advocated since several 
years by Neuberger and Narayanan at Nc —> oo, we are tempted to establish a 
connection with it : it could be a remnant of this at Nc = 3. We stress the fact 
that this property is revealed only through a polynomial extrapolation from non 
negligible quark masses of 30 MeV or more (bare masses). 
Alternative method for calculating the quark condensate As a by-product, our 
calculations suggest that the rather direct method of calculating the condensate 
through the chiral extrapolation of m J d A x < P^(x)P^(0) > is efficient. For our 
larger lattice, the results is very encouraging, considering the simplicity of the 
method : < ipip > I7g - (2 GeV) = -(2.2 ± 0.4) 10~ 2 GeV 3 . 
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A High orders of the Wilson coefficient of the 
condensate 

From the calculations of Chetyrkin and Maier ( [5] [5] ), we extract the 
following result at Nf — for Euclidean p : 

m(p 2 ) = 16w 2 /p 2 [a fl /7r(--L) + 
KAr) 2 (-|) + 

. , , 4 . 26331733 2236285 

(a s /n) 4 ( h C 3 

v ' ' v 497664 995328 sv ' 

-J^> + fli« 5 »i< 1 + s <' , » < 41 > 

a s being taken at p. It is easily seen that the first order is as given by Politzer, 
and the second one as given by Pascual and de Rafael. We have obtained these 
numbers by the following manipulations. The authors consider the renormalised 
OPE of the quark propagator renormalised in the MS scheme, and write the 
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contribution of the quark condensate as C(q) ^ < ipip >. They find for the 
scalar part and vector part of the renormalised propagator S — A{p) + B{p)ft : 

A(p) = -m/p 2 + C^{p) < ^ > +... 

B(p) = -1/p 2 + ... (42) 

where dots denote n > orders in a or power corrections with power larger than 
1/p 2 , or terms suppressed by n > powers of m. In A(p), the tree order term 
m is kept in a first step to indicate the conventions of the authors. We note 
that m{p 2 ) as we calculate from the ratio of bare quantities is also the ratio of 
the renormalised quantities in any scheme, since the factor Z2 applies to both 
terms of the fraction. Then, at m = 0, : 

m{p 2 ) = A(p)/B(p) = -p 2 C^(p) < ^ > (1 + E(p)) (43) 

where S(p) is the A/5 quark self-energy. Now, we pass to the Euclidean space 
by setting p 2 = —p 2 E , or making the substitution p 2 — ¥ —p 2 . Then C being read 
as l/p 4 C", where C" represents the radiative corrections and is a polynomial in 
a s , one ends on : 

m(p 2 ) = 1/p 2 C'{a s ) < W > (1 + E(p)) (44) 

C and < ■ipip > being negative, m(p 2 ) > as observed on the lattice. Finally, 
we set fi = p. Then E(p) is found to give a very small correction. We borrow 
the expression from Chetyrkin and Retey : [26] : 

1 + S(p) - 1 + l/4 2 (C A C F (41/4 - 3 C(3)) + C% (-5/8))(a s /7r) 2 
+l/4 3 (ClC F (159257/648 - 3139/24 C(3) - 69/16 C(4) + 165/4 C(5)) + 

CUCK-997/24 + 44 C(3) + 6 C(4) - 20 C(5)) + Cf,(-73/12))(a s /7r) 3 (45) 

or numerically : 

1 + E(p) = 1 + 1.59151(a s /7r) 2 + 23.2809(a s /7r) 3 (46) 
This is 1.005 for a s = 0.14 (p = 5 GeV) 
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